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Abstract 

It is shown that the scaling operators in the conformal limit of a two-dimensional 
field theory have massive form factors which obey a simple factorisation property 
in rapidity space. This has been used to identify such operators within the form 
factor bootstrap approach. A sum rule which yields the scaling dimension of such 
operators is also derived. 



1. The solution of conformal field theories represented a crucial step in our understanding 
of two-dimensional quantum field theories 0,0- It not only amounts to a complete 
description of the fixed points of the renormalisation group (RG), but also provides the 
starting point for the study of physical systems away from criticality. For example, one 
can consider the theory defined by the action 

A = Acft + 9 J d 2 x(p(x) , (1) 

as describing the perturbation of a conformal invariant theory by the operator (f(x) of 
scaling dimension 2 A < 2. Since the coupling constant g has physical dimension m 2 ~ 2A , 
m being a mass scale, the theory (Jl|) is no longer scale invariant; rather, it is associated to 
a RG trajectory flowing out of the original fixed point. In many cases an infinite number 
of integrals of motion survive in the perturbed theory and the resulting off-critical model 
is said to be integrable. A bootstrap procedure can then be applied, usually resulting 
in the determination of the exact particle spectrum and S-matrix of the theory || |4], |J. 
As a consequence of integrability, the S'-matrix turns out to be completely elastic and 
factorised. On the other hand, it is commonly believed that the knowledge of the S'- 
matrix amounts to a complete solution of a quantum field theory. In particular, it should 
encode the information about the operator content of the model and should enable the 
computation of correlation functions. The method which has proved so far very effective 
in dealing with such an ambitious program is known as the form factor bootstrap. Form 
factors (FF) are matrix elements of local operators 0(x) between asymptotic multiparticle 
states and will be denoted as [] 

Ff(e 1 ,...,e n ) = (o\o{o)\A(e 1 ),...,A(e n )) . (2) 

They are an interesting subject for theoretical investigation because their structure in- 
volves the particle description of the theory, an expression of the infrared dynamics, as 
well as the operator content, which is deeply related to the conformal structure of the ul- 
traviolet fixed point. Moreover, if the FF are known, correlation functions can be written 
down in terms of the spectral sum 

<0i(*)0 2 (o)> = £ ^ / ^ • • • f^?('x, • • • > <>») • • • . 0*)] * . (3) 

P£ being the total energy-momentum of the n-particle intermediate state. 

The determination of FF in integrable models proceeds through two basic steps. One 
first enforces the general constraints of analyticity, unitarity and crossing symmetry de- 
riving from the standard S- matrix theory || [7J. Since no information about the specific 

1 In order to avoid inessential complication of the notation, we refer throughout this letter to a theory 
whose spectrum consists of a single species A of particles with mass m. We also adopt the standard param- 
eterisation of the on mass shell momenta in terms of the rapidity variables Of. = (m coshOi, m sinhdi) . 
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nature of the operator 0(x) is provided at this stage, the general solution of the resulting 
system of functional and residue equations (linear in the operator) must correspond to a 
complete description of the operator content of the theory. It has been shown for several 
models that the dimensionality of the linear space of solutions of these equations coincides 
with that predicted by the underlying conformal theory [|3[, |9[, [T0[|. 

The second step consists in selecting out of the general solution the FF corresponding 
to particular operators, the scaling operators being the objects of main physical interest. 
Such operators provide the natural physical basis for the space of local operators of the 
theory. Nevertheless, it is by no means obvious how to identify this basis among the 
solutions of the FF equations. A first selection rule comes from the general constraint |T2 



lim F* (0i, . . . , 9 n ) < const. e A ^ e ' 1 , (4) 

relating the asymptotic behaviour of the matrix elements of a scaling operator to 
its scaling dimension 2A<j, Q. When the operator space splits into different sectors distin- 
guished by some internal symmetry of the theory, the asymptotic bound (01) is typically 
sufficient to fix the FF solutions. Yet, the problem remains conceptually interesting in 
absence of internal symmetries. We now show that in the latter case the FF of the relevant 
(A$ < 1) scaling operators of the theory are characterised by the following asymptotic 
factorisation property 

lim FX l (e 1 +a,...,e r +a,e r+1 ,...,e r+l ) = ^F^(e 1 ,...,e r )F^e r+1 ,...,9 r+l ) . (5) 

Such a property had been noticed in the past to be fulfilled by some FF solutions in 



specific models [11 



Consider the case in which the action ([!]) describes a massive integrable model without 
internal symmetries and denote by S(9) the two-particle scattering amplitude, 9 being the 
rapidity difference between the colliding particles. It can be argued on general grounds 
that 

lim S(9) = 1 , (6) 

and that this result implies that the FF of a relevant scaling operator actually tend 
to a constant (with respect to 9{) in the limit (f|), so that the limit in eq. (j|) is surely 
well defined. Also, due to the absence of internal symmetries, the vev ($) should be 
nonvanishing and can be written as 

($) = v$ m 2A * , (7) 

v$> being a dimensionless constant. 



2 We will refer to scalar operators. 
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The basic point to be realised is that the limit in the r.h.s. of eq. is actually a 
massless limit into the ultraviolet conformal point. In general, such limit is obtained 
sending to zero the mass of each particle while sending to infinity its rapidity in order 
to keep the energy finite jl^, [T3||. The resulting theory consists of right and left-moving 



massless particles, Ar and Al, whose Zamolodchikov-Faddeev operators may be formally 
obtained from those of the original massive particle A as 

A R (9) = \im a ^ +O0 A(9 + a/2) , (8) 
A L (6) = \im a ^ +00 A(9-a/2) . 

The dispersion relations are 

p° = p 1 = — e e for right-movers , (9) 

p° = —p 1 = —e~ e for left-movers , (10) 

M = me Q//2 being a finite parameter. The scattering amplitudes characterising the in- 
teraction of the masseless particles are readily obtained taking the limit of the massive 
amplitude: Srr(9) = Sll{9) = S(9) and Srl(9) = \im a ^ +00 S(9 + a). In the present 
context the latter limit for Srl gives a constant phase whose rapidity independence en- 
sures the decoupling of the right and left sectors and the scale invariance of the massless 
theory 0. In absence of internal symmetries, eq. (|6|) implies Srl = 1- In summary, the 
massless scattering theory defined through the above limiting procedure provides a par- 
ticle description of the conformal point to which the original massive model flows in the 
ultraviolet limit. 



Massless form factors can be introduced along the same lines ||15|| . In order to prevent 
a trivial vanishing of the matrix elements (0) when m — > 0, it is convenient to refer to the 
rescaled operator 

Since the matrix elements of a scalar operator depend on rapidity differences only, we can 
write 

lim Q ^ +00 Ff + i{9\ + a, . . . ,9 r + a, 9 r+1 , . . . , 9 r+ i) = 
lim Q ^ +00 Fj*. ,(0i + a/2, ...,9 r + a/2, 9 r+1 - a/2, 9 r+l - a/2) = 

F r l (9i, . . . , 9 r \9 r+1 , . . . , 9 r+ i) , (12) 

3 Whcn the right-left amplitude is rapidity dependent the theory describes a massless flow between 
two fixed points. 
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where 



...,9 r \9[,..., 6[) = (0m0)\A R (6 1 ), . . . , A R (9 r ), A L (9'i), A L {6[)) , (13) 

denotes a massless form factor. On the other hand, as a consequence of the decoupling 
of the right and left sectors at the conformal point, the massless form factor of a scaling 
operator can be written in the factorised form 

. . . AK, . . . , 9[) = 72*(0 a , ...,9 r ) Cf(9[, (14) 

Consider now a form factor in the massive theory and take its massless limit in which all 
the particles become right-movers. Due to Lorentz invariance, the rapidity shifts needed 
for the limit can be completely rescaled out and we have 

F*{0 U ...,9 r ) = F r %(9 u ...,9 r \) = K(9 U ■ ■ ■ , Or) 4 ■ (15) 

A similar equation holds when all the particles become left-movers. Since ($) = t>$ along 
the whole flow, at the critical point we have ^ 

111 = 4 = v'J 2 . (16) 

Using eqs. (|14|), (|15|) and fll6f) , eq. ( |12|) can be rewritten as 

lim Ff +l (6x+a, ...,9 r +a, 9 r+1 , . . . , 9 r+l ) = —F*(9 1} ...,9 r ) F*(9 r+1} . . . , 9 r+l ) , (17) 

in terms of massive FF only. Going back from $ to $, eq. (§) follows. 

It is worth noticing that eq. (j^) can be used to determine the vev (<&) from the knowl- 
edge of the multiparticle FF Q. Taking into account that only the vev of the perturbing 
operator can be obtained by other means (thermodynamic Bethe ansatz), this must be 
regarded remarkable circumstance. 

The crucial conditions entering the above derivation of the factorisation property (||) 
are the nonvanishing of all the FF of the considered scaling operator (including the vev) 
and their constant asymptotic behaviour. These features, which are guaranteed in absence 
of internal symmetries, are often shared by some operators in theories with symmetry. The 
FF of these operators then factorise asymptotically according to eq. (|). As an example, 
we mention the exponential operators in lagrangian theories like the Sinh-Gordon model 



11 ] , or the linear combinations a ± [i of the order and disorder parameters in the thermal 



4 Remember that $ is a rescaled (dimensionless) operator; the vev (Q) of the operator $ vanishes in 
the massless limit, as expected in conformal field theory. 

5 More precisely, what one can fix is the ratio which, being independent from the normalisation of 

the operator, is an universal number. 
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Ising model |TB| . It seems reasonable to expect that a suitable generalisation of eq. (|5|) 
should apply to any scaling operator even in presence of internal symmetries, but no 
general pattern has been identified so far. 



2. Once the FF solutions for the relevant scaling operators have been selected using 
the factorisation property (|5|), it remains to be established which scaling operator each 
solution corresponds to. For this purpose, it is important to be able to recover the 
scaling dimension 2A$ from the FF solution. In principle this can be extracted evaluating 
the short distance behaviour of the spectral series (||]) for the correlator ($(x)$(0)). In 
practice, better quantitative results can be obtained exploiting the properties of the stress 
tensor as the generator of dilatations. In order to illustrate this point it is useful to start 
with some perturbative consideration. 

The operator space of the conformal point and that of the perturbed theory have 
the same basic structure. In particular, a scaling operator $(a;) in the off-critical theory 
can be associated to a conformal operator Q(x) of scaling dimension 2A.J,. When doing 
that in a perturbative framework, however, renormalisation effects induced by ultraviolet 
divergences must be taken into account (nj. In fact, denote by X a generic product of 



operators and consider the usual perturbative expansion of the correlator 

(X $(0)) = (X $(0))cft + 9 [ d 2 x (X <£(0) <p(x)) C FT + 0(g 2 ) , (18) 

Je<\x\<R 

where the correlators in the right hand side are computed in the conformal theory, and 
e and R regularise the ultraviolet and infrared divergences, respectively. The integral in 
eq. fll8|) is UV divergent only if the conformal OPE 

<p{x) $(0) = £ CJ* |af (a*-a,-a) A k (0) (19) 

k 

contains operators with scaling dimension 2A& such that 

7 fc = A fc - A$ - A + 1 < . (20) 

In this case we obtain a first order UV finite correlator (X$(0)) by defining the renor- 
malised operator as [] 

<S> = $ + gJ2he 2 ~< k A k + 0(g 2 ) , (21) 



where 



(22) 



Ik 



3 Eq. (EH]) must be suitably modified in the case of logarithmic divergences (7^ = 0). 
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and the sum runs only over the operators satisfying the condition (^). This implies 
that, in general, renormalisation mixes the original operator $ with a finite number of 
operators of less scaling dimension. 

Consider now the euclidean correlators of the operator $ with the components of the 
stress energy tensor T = \{T\\ — T22 — 2zTi 2 ) and 9 = T\\ + T22 

Fizz) 



<T(M)*(0)> 



z 
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(6(^)$(0)> C =^M, (23) 

where z and z are complex coordinates. Conservation of the stress energy tensor espressed 
by the equation 

dT + ^dO = , (24) 

leads to the differential relation 

D = \G , (25) 

where D = F + jG and the dot stays for zz-A=. Since the trace of the stress tensor is 
related to the perturbing field as 

e(x)=47r<7(l-A)<^r) , (26) 

the short distance behaviour of the function G is determined by the OPE (W) 

G{x) ~ 2^(2 - 2A)C ^ \x\ 2 ^ (Ac) , x - (27) 

where we denoted by A the most relevant operator appearing in (|T9|). We can now 
distinguish two basic cases: 

a>) 7o > 0. In this case G vanishes as x —>■ (conformal limit) and we conclude that 
the function D is stationary and coincides with F at the fixed point. Since the operator 
does not mix under renormalisation, its conformal OPE with T(x) can be safely 
used in a neighbourhood of the fixed point 

F(x) ~ A $ ($) , x -> . (28) 

If the theory described by the action (|l]) corresponds to a massless flow between two fixed 
points, a similar analysis can be repeated in the neighbourhood of the infrared fixed point. 
Then, integrating eq. ( p5| ) over all distance scales, one finds 

- Ai R = J d 2 x (0(x) $(0)) c . (29) 

In a massive theory A$ R = 0. We will illustrate in a moment with few examples the 
effectiveness of this sum rule within the FF approach. 
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b) 7o < 0. The function G no longer vanishes at x = and an attempt to use the sum 
rule (|29|) would be frustrated by the divergence of the integral. This is a consequence of 
the fact that the operator &(x) now mixes under renormalisation (see eq. ([H])) so that 
the function F no longer behaves as in ( pBD at short distances. The correct behaviour can 
instead be obtained integrating eq. (|25|) . One finds 

F(a;)^7r^(l-A)C^^— ^|a;|*»(i4o> , 7o < (30) 
* 7o 

2^(1 -A)Cj $ log |x|(A ) , 7o = 0. (31) 

Actually, eq. (|29|) is nothing but the Ward identity expressing the fact that is the 
generator of scale transformations, and as such a suitable generalisation is valid in arbi- 
trary dimension. This may be derived simply by modifying the arguments presented, for 
example, in ref. Consider the effect of making a nonuniform infinitesimal RG trans- 
formation, which generates the coordinate change x M — > x ni = x^ + where 5x^ = ex M , 
corresponding to uniform scale transformation, in the region between two hyperspheres, 
of radii R\ and R2, while 5x^ = in the regions \x\ < R[ < R\ and \x\ > R' 2 > Ri- 
Within the thin shells R[ < \x\ < R\ and R2 < \x\ < R' 2 , Sx^ is chosen to be an arbi- 
trary differentiable function which matches smoothly onto the other regions. From the 
definition of the stress tensor, the change in the hamiltonian (action) is 

6S = -(1/S d ) J VSaTT^dk, (32) 

where Sd, the area of a d- dimensional sphere, is conventionally included. This gives a 
contribution /jj 1 <| a; | < ij 2 OgHe from this region. Within the two shells, integrate by parts: 
the bulk term vanishes by the conservation d^T^ = of the stress tensor, leaving two 
surface integrals on \x\ = R\ and R2 of the form / x^T^ v dS u ' . Now consider evaluating 
($(0)) with respect to this modified hamiltonian. Since no change has been made for 
\x\ < R\, this quantity is in fact invariant. Thus the first order changes from the bulk 
and boundary terms should cancel. Expanding these out from the hamiltonian to O(e), 
the bulk term gives 

-4-/ <e(*)$(0)) c ^, (33) 

&d JRi<\x\<R 2 

while the boundary terms involve integrals over (T rr .(x)$(0)). In the limits -Ri — > and 
i?2 —> 00 these may be evaluated using the operator product expansion formulae pO 
appropriate to the UV or the IR conformal theories respectively: 

T^MO) = te*x^ v ~^lf x2 ^ m + . . . , (34) 
where x$ is the scaling dimension of $. This then leads to the result 

(^-4^(4) = "I" / (e(x)$(0)) c d*r, (35) 
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which is the generalisation of (P§[), with x<j> twice the complex dimension A$. 

In the case when the right hand side diverges, it is not permissible to take the limit 
— » 0, i?2 — ► oo), and, as above, this may be seen through the operator mixing which 
is exhibited in the operator product expansion fl3~j|). 



3. The Ising model provides an interesting example for testing the ideas discussed in this 
letter. The situation is particularly simple in the purely thermal case, which in terms 
of the action (P corresponds to perturbing the Ising critical point by the energy density 
operator e(x) of scaling dimension 2A £ = 1 (g ~ T — T c ). The only other relevant operator 
in the model is the local magnetisation er(x) of scaling dimension 2A CT = 1/8. The theory 
can be described in terms of a free majorana fermion of mass m, corresponding to a 
scattering amplitude S(6) = —1. The only nonvanishing FF of the components of the 
stress tensor T and G = 2irge (bilinear in the fermions) are 

6>i — 9 2 



iff 



h,o 2 ) 

h,9 2 ) 



-2mm 2 sinh ■ 



m 2 e 8l+62 sinh 



2 
Ox 



9, 



2 2 

The correlators of a with T and vanish for symmetry reasons at T > T c . Alternatively, 
one can work at T < T c , or exploit duality and refer instead to the disorder operator \i. 
Here we only need its two-particle FF 

6*i — 6 2 



#2) = tanh ■ 



(36) 



The following (euclidean) correlators are then obtained using eq. 



(T(x)T(O)) 

(T(x)e(o)) 

(T(*)M0)> 

<e(*)e(o)) c 
(e(x)/i(o)) c 



m 

T6 

4 - 
m z 

(A*> 



i^i(m|x|)-K"3(m|x|) — i^2( m l a: l) 
^(m|x|) — K (m\x\) K 2 (m\x\) 



16z 2 



-2m\x\ 



m 



K 2 {m\x\) — KQ(m\x\) 

^—2m\x\ 

— - + Ei(-2mb|) 

2mx v 1 w 



-m (/i) 



where K v are the modified Bessel functions and Ez the exponential integral function 
Ez(— x) = — dt exp(-t) jt . One can easily check that the above correlators satisfy the 
differential equation (p5|). Comparison with the short distance expectations 



(T(x)T(O)) 



c 

2z* 
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(T(x)e(O)) ~ ^(l-A £ )log|x| , 

(T{x)n(0)) ~ , \x\^0 

z z 

(e(x)e(O)) ~ |x|- 4Ae , 
(e(xMO)) * C £ »N- 2A % 

gives c = 1/2, A M = 1/16, A e = 1/2, m = 2irg and = —1/2. In writing the second 
equation above we took into account that 7o = for e and used eq. (^) with A = I. Of 
course A M can also be obtained inserting the exact correlator (G(x)^(0)) c into the sum 
rule (H). 

Less simple is the case of the purely magnetic perturbation of the Ising critical point 
{<p(x) = o~(x) i n the action (|l|)). The theory is again integrable but the spectrum now 
consists of eight massive particles Q. Since the presence of the magnetic field breaks 
the invariance under spin reversal, no internal symmetries are left in the model. The FF 
bootstrap program described in the first part of this letter was carried out for this model 



in ref. [p~2"| , [HJ . In particular, the factorisation equation (|5]) was used in ref. |21| to identify 
the FF solutions for the two relevant scaling operators a and e. Their scaling dimensions 
were also effectively estimated using the sum rule (|29|). Applications of the factorisation 
property @ to other models are presented in ref. . 

Finally, as an example of the massless case, let us consider the flow from the tricritical 
to the critical Ising points (the minimal models M.^ and of conformal field theory, 
respectively). The model is integrable and can be described in terms of a single species of 
massless particles whose scattering is characterised by the amplitudes Srr = Sll = — 1 
and Srl(6) = tanh (f - f ) JEJ. The correlator (6(ar)^(0)) c can be computed through 
the spectral series (|3p using the results of ref. Jl5|, where FF of several operators of the 
theory were determined. Here we again consider for symmetry reasons the disorder pa- 
rameter n instead of the magnetisation a. Both operators have scaling dimensions 3/40 in 
the ultraviolet limit and 1/8 in the infrared limit. Being the most relevant operator in the 
theory, /1 does not mix under renormalisation and the sum rule ( p9f) can be safely used to 
evaluate the total variation in its scaling dimension along the flow. The integration of the 
first (four-particle) contribution gives —0.0255; the addition of the second (six-particle) 
contribution leads to the result —0.0249, showing that the FF series rapidly converges to 
the expected result A^ v — A IR = 3/80 — 1/16 = —0.025 in spite of the massless nature 
of the theory. 
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